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ABSTRACT 



Aims. We study the variation in the magnetic field strength and the umbral intensity of sunspots during the declining phase of the 
solar cycle no. 23 and in the beginning of cycle no. 24. 

Methods. We analyze a sample of 183 sunspots observed from 1999 until 2011 with the Tenerife Infrared Polarimeter (TIP) at the 
German Vacuum Tower Telescope (VTT). The magnetic field strength is derived from the Zeeman splitting of the Stokes- V signal 
in one near-infrared spectral line, either Fe 1 1564.8 nm, Fei 1089.6 nm, or Si 1 1082.7 nm. This avoids the effects of the unpolarized 
stray light from the field-free quiet Sun surroundings that can affect the splitting seen in Stokes-/ in the umbra. The minimum umbral 
continuum intensity and umbral area are also measured. 

Results. We find that there is a systematic trend for sunspots in the late stage of the solar cycle no. 23 to be weaker, i.e., to have a 
smaller maximum magnetic field strength than those at the start of the cycle. The decrease in the field strength with time of about 
94Gyr~' is well beyond the statistical fluctuations that would be expected because of the larger number of sunspots close to cycle 
maximum (14 Gyr -1 ). In the same time interval, the continuum intensity of the umbra increases with a rate of 1.3 (± 0.4)% of I c yr" 1 , 
while the umbral area does not show any trend above the statistical variance. Sunspots in the new cycle no. 24 show higher field 
strengths and lower continuum intensities than those at the end of cycle no. 23, interrupting the trend. 

Conclusions. Sunspots have an intrinsically weaker field strength and brighter umbrae at the late stages of solar cycles compared 
to their initial stages, without any significant change in their area. The abrupt increase in field strength in sunspots of the new cycle 
suggests that the cyclic variations are dominating over any long-term trend that continues across cycles. We find a slight decrease in 
field strength and an increase in intensity as a long-term trend across the cycles. 

Key words. Sun: sunspot - Sun: magnetic fields - Sun: photosphere - Sun: evolution 



1. Introduction 

Morphological and physical studies of sunspots have b een de- 
scribed for several decades dBrav & Loughheadlll964l) . It was 
already known half a century ago that lar ger umbrae are darker 
and sh ow higher magnetic field strength. lAlbregtsen & Maltbvl 
d!978l) first reported a systematic variation in the umbral in- 
tensity in the infrared as a function of the solar cycle, with 
dark er umbrae appearing in the early phase of the cycle (see 
also lAlbregtsen & M altbv 198lh. lAdiabshirzadeh & Koutchmvl 
(Il983l) attributed this variation to the change in the number 
of umbral dots within the solar cycle. This was dismissed by 
lAlbregtsen et al.l (|1984|) . who showed that the measured varia- 
tions cannot be solely caused by umbral dots. Subsequent stud- 
ies of the umbral intensity and the magnetic field strength es- 
tablished a relation between the continu um intensity and the 
mag netic field strength in sunsp ot umbrae dKopp & Rabinll992t 
Mar tinez Pillet & Vazquezl[l993l) : the greater the field strength, 
the lower the intensity. While new high-resolution observations 
are still consistent with this empirical relation, no clear exp lana- 
tion has been provided (see, however. ISchtissler et al.lll994h . 

The variation in the field strength of sunspots in the 
last two cycle s has been the subject of several stud- 



ies (ILivingstonl . 

I2004t iPenn & MacDonald l2007t iLeonard & Choudhary 2008; 



2002; Pennetal. 2003; Norton & Gilman 



ISchad & Pennl l2010l) . ILivingstonl d2002l) measured the sepa- 
ration of the cr-components of Stokes-/ of the Fei line at 
1564.8 nm, which is a Zeeman triplet with an effective Lande- 
factor of 3.0, and concluded that sunspots in cycle no. 23 were 
intrinsically weaker than those in cycle no. 22. A similar data 
set was also used in subs equent works dLivingston et al.l 120061; 
IPenn & Livingston] 1201 lh . These authors found that the field 
strength of sunspots was stronger at the beginning of cycle no. 22 
than at its end. Although th e magnetic sensitiv ity of the infrared 
Fe 1 1564.8 nm line is high (iRtiedi et alJI 19951) . the influences of 
scattered light and line blends in the intensity profiles on the 
measurement motivates one to employ more accurate methods 
to derive the field strength. It is not clea r how dynamo the ory 
can explain such cyclic dependency (e.g.. lOssendrrjveil2003l) . 

In this contribution, we used spectropolarimetric data in the 
Fe 1 1564.8 nm line and other near-IR lines to evaluate the cyclic 
variation in the magnetic field strength of sunspots. We used full 
Stokes spectra to determine the magnetic field strength rather 
than solely using the intensity profiles. The advantage with re- 
spect to the stray light contamination is that the intensity in 
the quiet Sun (QS) surroundings is higher than in the umbra, 
whereas exactly the opposite relation holds for the polarization 
signal. Therefore intensity spectra in the umbra are contaminated 
by contributions from the QS, but the splitting of Stokes- V is not. 
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Table 1. Atomic 
lines dNave et al.lfl99 



roperties of the observed spectral 



Line A (nm) Exc pot (eV) log(g/) g-effective 



Fei 
Fei 

Si i 



1564.852 
1089.630 
1082.709 



5.426 
3.071 
4.954 



-0.669 

-2.845 
0.363 



3.00 
1.50 
1.50 



Another advantage of this sample compared to the previously 
mentioned works is that the majority of sunspots were observed 
with an active image correction; i.e., the spatial resolution of our 
data set is better than observations prior to the advent of correla- 
tion trackers or adaptive optics (AO). 

Section|2]describes the observations used. The data analysis 
methods are explained in Sect. [3] Section|4]gives the results, to- 
gether with numerical tests. The findings are discussed in Sect. [5] 
Section[6]provides the conclusions. 



2. Observations 

We analyzed a sample of 231 sunspots observed from 1999 to 
2011 with th e Tenerife Infrared Polarimeter (TIP-I and TIP- II , 
respectively; iMartinez Pillet et al.l 119991: ICollados et"ai] 120071) . 
All observation were carried out at th e German Vacuum Tower 
Telescope (VTT. ISchroter et aDll985l) . The observed maps usu- 
ally cover a sunspot and the surrounding QS area. For some of 
the data taken with TIP-I using its small field of view (FOV) 
of about 30", the sunspots' penumbra was not fully covered. 
Each data set consists of full Stokes profiles in magnetically 
sensitive infrared lines such as Fei 1564.8m, Fei 1089.6 nm, or 
Si 1 1082.7 nm. TableQ]summarizes the atomic parameters of the 
selected spectral lines. Out of the 231 maps, 99 were observed 
at 1.56yum and 84 at 1.1 fim. The remaining 48 maps mainly 
come from the earliest observations with TIP-I and were taken 
in some uncommon wavelength ranges, often covering molecu- 
lar lines. These data could therefore not be used without consid- 
erable effort. We thus only selected those observations in which 
one of the three spectral lines listed in Table[T]was recorded. That 
amounted to 183 full Stokes sunspot maps covering the descend- 
ing phase of cycle no. 23 and the rise of cycle no. 24. 

The low temperature of umbrae allows certain molecules to 
form. As a result, the umbral profiles are usually contaminated 
with molecular blen ds (e.g., the FTS atlas of umbral profiles of 
I Wallace et al.l20 00ii). FigureQ]shows examples of umbral Stokes- 
/ and V profiles at 1.56 /mi for two different heliocentric angles 
(fi - 0.58 and 0.76). The difference in the continuum level of the 
Stokes-/ profiles is real. As seen in this figure, there are sev- 
eral strong blends around the Fei 1565.2 nm line that corrupt 
not only the intensity profile but also the Stokes- Q/U/V pro- 
files of this line (bottom panel, Fig. Q]). For this reason, we de- 
cided not to use the Fei 1565.2 nm line. The molecular blends 
are, however, useful for determining of a tmospheric strati fica - 
tions in inversions of umbral spectra Ce.g.. lMathew et aLll2003h . 
Ther e are several blends aro und the Fe 1 1564.8 nm line as well 
(e.g., iLivingston et al.ll2006l their Table 1). The main blend in 
this line is a CO line to the blue of the line center wavelength. It 
is, however, weaker than the OH blends near the Fe 1 1565.2 nm 
line. 

The spatial and spectral resolution, as well as the rms noise 
level, is different for each map. A careful analysis of the individ- 
ual sunspots was done to take the diversity of the data into ac- 
count. While the data until 2005 were recorded using a correla- 
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Fig. 1. Sample of two umbral profiles. The solid (dashed) line 
corresponds to a sunspot at fi — 0.58 (0.76). The continuum in- 
tensity of the quiet Sun profiles corresponding to each of them 
is unity. The short vertical lines mark four OH spectral lines, 
while the long vertical line denotes the Fe 1 1564.8 nm line. Top: 
Stokes-/, bottom: Stokes- V. 



tion tr acker for real-time image correction (CT, iBallesteros et al.l 
1996), the Kiepenheuer-Institute adaptive optics system (KAOS, 
Ivon der Liihe et all .2003) was utilized for all later data to im- 
prove the spatial resolution. 

The old TIP-I camera dMartinez Pillet et al.|[l999h was used 
until mid 2006, providing a slit length of 27 arcsec at a spatial 
samplin g of 0736 per pixel. The integration of the new TIP-II 
camera (ICollados et al.ll2.007h increased the slit length and at the 
same time doubled the spatial and spectral resolution because 
of the larger CCD chip and the smaller pixel size. Examples of 
data sets with the old and new TIP systems are shown in Fig. 
|2] The top panel shows a TIP-I data set with CT correction, 
while the bottom panel shows a data set recorded with TIP-II 
taken with AO correction. An overview of the data up to 2009 
is available online in the TIP archivqj. The header of each file 
contains information about, among other things, the location of 
the spot, scanning steps, and the integration time. This allowed 
us to retrieve the position and dimension of the sunspots. We 
used this information to calculate the area of the umbra and the 
sunspot as a whole. The calibration of the data as well as the 
cross-talk correction was perform ed in the same manner as in 
ISchlichenmaier & Colladosl d2002l) . 



1 http://www3.kis.uni-freiburg.de/~cbeck/TIP_archive/ 
TIP_archivemain.html. If you are interested in using the data, please 
contact mcv@iac.es. 
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Fig. 2. Example of TIP scans. Top: sunspot on 24 May 2001, 
recorded in Fe 1 1089.6 nm with CT correction. Bottom: sunspot 
on 17 July 201 1, recorded in Si 1 1082.7 nm with AO image sta- 
bilization. Tick marks are in arcseconds. 

3. Data analysis 

For each data set, an average quiet Sun profile was calculated 
using a region in the FOV outside of the penumbra. We used 
two of the spectral lines covered in each wavelength range to 
determine the spectral dispersion for each map. Then, the rms 
noise of the polarization signal, cr, was estimated in a contin- 
uum band. After that, the amount of separation between the 
two lobes of Stokes-V was measured for all locations on the 
map with a maximum Stokes- V amplitude larger than 12 cr. 
"Lobes" were defined as local maxima in the polarization sig- 
nal of Stokes-V that a) were larger than the five wavelength 
pixels to the left and right and b) had an ampli tude above the 
12-cr level (cf. lBeckll2006t iRezaei et al.ll2008h . This conser- 
vative threshold ensures that we only analyze clear antisym- 
metric Stokes-V profiles with a drawback that we have some 
gaps in the maps in the neutral lines of Stokes-V where the po- 
larization signal strongly decreases. Because there are only a 
few data sets close to the limb (p <0.5), the neutral line only 
passed through the umbra in a couple of cases. In particular, 
the missing values in the neutral line do not affect our results 
concerning the maximum magnetic field strength in umbrae. 
Whenever the profile exhibited more than two lobes, the two 
lobes most separated in wavelength were selected. Finally, we 
constructed maps of the magnetic field strength from the deter- 
mined location of the two lobes in the Stokes-V profi les and the 
spectral dispersion in the strong field a pproximation (Stix 2002; 
iLandi degl'Innocenti & L andolfi 2004): 



B = 



A,! 



4.67xl0- 12 ^ 2 geff 



(1) 



where A and A/1 are in nm and B is in Gauss. In the strong field 
limit, the separation of the two lobes is independent of the incli- 
nation angle, hence of the location of the spot on the disk. The 
signal amplitude is, however, affected by the inclination angle, 
but we do not use it here. In a plot of the maximum field strength 
v*. cosine of the heliocentric angle (not shown here), we did not 
find any significant trend. 
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Fig. 3. Variation in the minimum continuum intensity in the 
umbra as a function of maximum magnetic field strength. The 
1.1 /mi data was multiplied by a factor 1.3. Red and black sym- 
bols mark the 1.1 and 1.6jum data, respectively. P-incomplete 
and U-incomplete mark sunspots whose penumbra and umbra 
were not fully covered in the map. The solid line shows a 
least-squared fit, and the dashed lines indicate the \cr interval. 
Numbers in parenthesis denote the 1 cr uncertainty. 



The uncertainty of the estimated field strength was calcu- 
lated using standard error propagation applied to Eq. (|T). The 
main source of uncertainty stems from the measurements of 
the separation of the lobes, whereas the effective Lande fac- 
tor and the rest wavelength are known with high accuracy. We 
assumed an error of 0.1 pixels in the derivation of the posi- 
tion. This translates then into a (random error) ctb of about 
20 to 60 G, depending on the wavelength range and the disper- 
sion of the spectra. The value is comparable to the formal er- 
ror attribut ed to field strengt h in inversions of the Fe i line at 
1564.8 nm dBecket al.ll2007l) . There are also systematic errors 
caused by gradients of the field strength with height in the atmo- 
sphere. With a typical measured vert ical gradient of the magnetic 
field strength of 1 - 2Gkm" 1 (e.g. JSchmidt & Baithasailll994t 
IWestendorp Plaza et al.ll2001T) . this systematic error amounts to 
some 20 - 40 G difference in the measurements between p = 1 
andyU = 0.5. 

Besides the maps of the magnetic field strength, we also de- 
termined the continuum intensity and the area of each sunspot. 
To define the umbra and penumbra, we produced at first an 
umbral mask using a variable threshold in continuum inten- 
sity. Using a fixed threshold did not produce satisfactory re- 
sults because the umbral contrast varies w i th wavelength (e.g., 
Mattig j ll97U lAlbregtsen&Maltbvlll98U IChap man & Meverl 



1981 



Maltbvet alj 119861: iTritschler & Schmidtf 120021) . Using 
those (initial) masks, we constructed contours for the umbra and 
penumbra manually. Light bridges were excluded from the um- 
bral masks. In some data sets, such as shown in the top panel of 
Fig. |2 the penumbra was not covered completely, so the penum- 
bral area and the total sunspot area are underestimated in a few 
cases. The continuum intensities were not corrected for the limb- 
darkening effect. There was no significant trend in the minimum 
continuum intensity vs. the cosine of the heliocentric angle. 

A code was attributed to each spot describing its type in 
terms of complexity, completeness, and similar morphological 
characteristics (see Sect. 14. It . In the case of complex active re- 
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Fig. 4. Variation in the minimum continuum intensity in the um- 
bra as a function of area. Top: total sunspot area, bottom: umbral 
area. The solid lines indicate a bisector least-squared fit to all 
data points expect complex sunspots (circles). The dashed lines 
show the lcr interval. Symbols are like Fig. [3] 



gions (e.g., the top panel of Fig. |2), only the largest spot in the 
map was selected. 

The continuum intensity sample is not fully homogeneous 
because we have the continuum intensity in two different wave- 
length windows (1.1 and 1.6 /mi). As mentioned above, sunspots 
have higher contrasts at shorter wavelengths. We multiplied the 
continuum intensity of sunspots at 1 . 1 fim with a factor of 1 .3 to 
have the same scale as in 1.6yum, both in scatter plots and his- 
tograms (Sect. @}. This coefficient of relative umbral intensities 
of 1.3 corresponds to an umbral temperature of about 4200 K. 
The intensities in each wavelength range were normalized to the 
intensity in the QS surroundings of the sunspot. For a tempera- 
ture of about 6000 K in the QS, the ratio of continuum intensities 
in the Planck curve is /qs(/1 = 1100nm)//Qs(/l = 1565 nm) ~ 
2.77; i.e., the intensity at the shorter wavelength has to be mul- 
tiplied by 2.77 to be normalized to Iqs(A = 1565 nm). This 
yields a relative umbral intensity I um bm(A - 1100nm)//Qs(/l = 
1565 nm) = 1/1.3 • 2.77 ~ 2.13 , when the umbral intensity 
at 1565 nm is set to unity. A ratio of 2.13 in the intensities 
I(A = 1 100 nm)//(i = 1565 nm) is obtained in the Planck curve 
for T — 4165 K. The value fits to previous determinations of um- 
bral temperatures (e g.,[M altbv et al. 1986; Collados et al 1 1994: 
Mathew e t alJl2003l:ISolankil2003l:ISanchez Cuberes et al.l2005t 
Beck 2008). 



4DDD 


. . . 


i 


■ ^ ^ ■ j^- ■ ■ ...- 






-slope 


= 0.93 G arcsec" 2 


V >^ r" 


■ 


35DD 




-'a 

3 D D v<^v* 


?'""' 


- 


g" 3DDD 

X 

Hi 

E 
rn 25DD 


Co 

-a 






~ 


2DDD 


-jy' 






- 


15DD 




i .... i .... i 




■ 



5DD 1DDD 15DD 2DDD 25DD 

Sunspot area |aicsec 2 | 



4UUU 


-slope =3.32 G arcsec 2 






^ 


35DD 




-^ 


jr^"^ a j-" 


'- 




o ^ ■* ^^ 






. 


g- 3DDD 


to i°^i aj^ 






- 


115 

E 








■ 


03 25DD 


vv$^£fof& Di o 




Q 


- 


2DDD 


v v^i^ g a 




o 


: 


15DD 


■ i 


. . . i . 




■ 



1DQ 2DD 3DD 4DD 

Umbral area [arcsec 2 ] 



5DD 



Fig. 5. Variation in the maximum magnetic field strength in the 
umbra as a function of area. Top: total sunspot area, bottom: um- 
bral area. Symbols are like Fig|3] Linear fits are similar to Fig. [4] 



4. Results 

4. 1 . Empirical relations 

Figure[3]shows the scatter plot of the minimum umbral intensity 
vs. the maximum field strength, B max . By "minimum" we mean 
the average of the ten pixels with the lowest intensity. Similarly 
by B max , we mean the average of the ten pixels with the largest 
field strength. These pixels partly overlap those with the lowest 
intensity. Different symbols in the plot indicate morphological 
details of each sunspot. Triangles denote sunspots whose penum- 
bra touches the map border, although they have a complete um- 
bra. Circles mark sunspots in complex active region. This group 
of sunspots sometimes have extended or incomplete penumbrae. 

As a result, a triangle, for instance, does not exclude the spot 
belonging to a complex active region. Reverse triangles mark 
those objects that are either a pore or more like a protospot: they 
do not have a fully fledged penumbra. In rare cases, the map 
does not fully cover the umbra. These objects are marked with 
asterisks. As seen in Fig. [3] there is a tight correlation between 
the continuum intensity and the field strength. The classical or 
Pearson correlation (C p ), as well as the nonpa rametric Spearma n 
correlation coefficients (C s ), were calculated dPress et al.lll992l) . 
We find correlation coefficients of C p = 0.80 and C s = 0.76. The 
solid line in Fig.[3]is a standard least-squared fit with an slope of 
-0.19 (+ 1.5 x 10~ 2 ) of I c kG l and an offset of l.ll(+0.04) of 
I c . The offset of 1 1 % (solid line, Fig. |3]l above unit continuum 
intensity for BeOG matches stray light estimates in the umbra 
well, where stray light levels between 8 % and 15 % were found 
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Fig. 6. Variation in the maximum magnetic field strength (left), minimum continuum intensity (midddle), and umbral area (right) vs. 
time. Symbols show the individual sunspots while the blue filled circles mark the annual averages. The linear fits was performed for 
individual data points between the two vertical lines. The dashed vertical lines mark maximum and minimum of cycle 23. Symbols 
are like in Fig [3] 



for si milar data (e.g., iRezaei et al.ll2007t lBeckll2008b iBeck et all 
l20ll . 

Figure [4] shows the variation in minimum continuum inten- 
sity vs. area. In this figure as well as Fig. [5] we do not show the 
error bars because they are usually smaller than the size of the 
symbols. There are a few points that are off the main trend. A 
close inspection showed that those sunspots belong to complex 
active regions and have extended penumbrae. If we omit their 
extended penumbrae, their areas shrink significantly. As seen in 
Fig. [4] there is an anticorrelation between the minimum contin- 
uum intensity and umbral (sunspot) area. The correlation coeffi- 
cients are C p = -0.43 and C s = -0.40 (C p = -0.48 and C s = -0.57), 
respectively. This is consistent with earlier findings: larger um- 
brae are intrinsically cooler. The linear fits in Figs. 5l an d|5]are 
bisector fits. That means we assume that both variables are inde- 
pendent dAkritas & Bershadvll996l) . For these four fits, the com- 
plex sunspots (circles) were excluded. A standard least-square 
fit results in shallower slopes. Assuming a linear relation be- 
tween the B max and minimum intensity with the sunspot/umbral 
area is equivalent to the relation used by other authors (e.g., 
ISchad & Pennl2010l) . who fit a quadratic function between these 
quantities and the umbral radius. 

Figure shows the variation in the B max vs. area. The trend 
is clearer with the sunspot area than the umbral area. We find 
a positive correlation between the maximum field strength and 
the umbral (sunspot) area (Fig. [5). ForB max vs. umbral (sunspot) 
area, we find C p = 0.26 and C s = 0.20 (C p = 0.25 and C s = 0.39), 
respectively. If we skip the complex sunspots with the largest 
area, the correlation coefficient increases significantly. We cre- 
ated two random vectors of the same size as our TIP sample 
and calculated their correlations. After repeating this experiment 
for 1000 times, the standard deviation both for the Pearson and 
Spearman correlations were ctq = 0.07, so the found correla- 
tion coefficients are well above the distribution of purely random 
numbers. 

Part of the scatter in Figs. [4]and[5]is real because there is a 
diversity in properties of simple and complex sunspots, as well 
as differences betwe en fully developed and forming sunspots 
dCollados et al.lll994b . As discussed in Sectj3] the separation of 
the two lobes of Stokes- V in the strong-field regime is not a func- 
tion of the inclination angle nor the heliocentric angle. 



4.2. Temporal evolution 

Figure [6] shows the temporal evolution of the maximum field 
strength, the minimum continuum intensity, and the umbral area. 



Symbols display individual umbrae while the blue circles mark 
annual averages. The two dashed vertical lines at 2001 and 2009 
denote the maximum and minimum of the solar cycle no. 23 
(defined from the sunspot number). A decrease in the magnetic 
field strength is seen in the diagram (left panel). We examined 
correlations to check for any systematic trend of B max with the 
solar cycle. To this extent, we selected the time span between 
the two dashed lines in Fig. [6] The correlation coefficient for the 
magnetic field strength is about C p = -0.57 and C s = -0.43 for the 
Pearson and Spearman definitions, respectively. This indicates 
the significance of the decrease in magnetic field strength in the 
declining phase of solar cycle no. 23. 

We used a Bayesian method (e.g. iGregorv & Loredolll992t 
lGregorv||2005l) to fit a straight line to the individual data points 
between 2001 and 2009. To this end, a Markov Chain Monte 
Carlo (MCMC) m ethod was utilized with one million iterations 
dPatilet alj|2010h . This resulted in a linear fit with an offset of 
2915( + 30)G and a slope of -94(±7)Gyr 1 . The solid line in 
Fig. [6] shows the best-fit line. The dashed line shows a standard 
least-square linear fit whose slope is -88 (+ 14)Gyr -1 . The ex- 
tension of the solid line (Fig.[6]l to 2012 demonstrates that the 
decrease in the field strength does not continue in cycle no. 24. 
If we use the range of 2000 to 2009 instead, the slopes drop 
to -59 ( + 8) Gyr 1 and -56 (+ 14)Gyr~' for the Bayesian and 
least-square methods, respectively. 

As seen in the middle panel of Fig. [6] there is a weak 
tendency in sunspots such that the umbra is brighter at the 
end of the cycle. Both linear fits have a common slope of 
1.3(±0.4)%I c yr _1 . From the relation between the umbral in- 
tensity and B max (Sect. 14. Il l, one expects to find an increase in 
the continuum intensity because we find a decrease in the field 
strength. The amount of decrease of B max during the declining 
phase is about 750 G, leading to a predicted increase in the mini- 
mum intensity of some 14 % of I c . As a result, the expected slope 
for the temporal evolution of minimum intensity based on the 
empirical relation between B max and intensity (Fig. [3} is about 
1.8 %yr -1 , consistent with the linear fit within the lcr uncer- 
tainty range. 

The correlation coefficient of the umbral area vs. time (right 
panel, Fig. [6) in the same time interval is C p = 0.05 and C s = 0. 12 
for the Pearson and Spearman correlation coefficients, respec- 
tively. It stays insignificant for the deprojected area of the umbra. 
The corresponding coefficients for the sunspot area are C p = - 
0.06 and C s = -0.03, respectively, so such correlations are not 
significant (recall the limit of <x c =0.07 for purely random vari- 
ables, Sect. 14. 11 1. That means that within the scatter and statistics 
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Fig. 7. Probability density function of the maximum magnetic 
field strength. The red curve is the best-fit normal distribution to 
the data using a Bayesian method. 
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Fig. 8. Yearly sunspot number. The solid line shows the SIDC 
annual sunspot number. The filled circles mark the number of 
TIP maps in each year. 



of our data, there is no systematic change in the umbral size with 
the solar cycle. 

4.3. Distribution of parameters 

The probability density function (PDF) of the maximum field 
strength is shown in Fig. [7] A Kolmogor ov-Smirnov test w as ap- 
plied to find the best PDF for the data dPress et al.lll992h . This 
test allows one to reject that the observed and fit distributions 
are significantly different. A normal distribution with an average 
field strength of 2.6 kG and a variance of 0.3 kG reproduces the 
observed PDF of the field strength well. There is no umbra with a 
magnetic field strength below about 1 .8 kG The assumed normal 
distribution thus exceeds the observed PDF for field strengths 
below that value. In the limit of extremely large field strengths 
above about 3 kG, the best-fit PDF also does not matches the ob- 
served one perfectly. This will presumably be caused by the poor 
statistics of those few umbrae with a maximum field strength 
above 3 kG in the observations, and the selection of the largest 
sunspot in complex active regions. 

4.4. Numerical experiment 

To investigate the influence of noise and systematic effects on 
the derived slope of B max with the phase of the solar cycle, we 
performed three numerical tests. This allowed us to evaluate the 
amount and importance of spurious signals in our results. In par- 
ticular, we tried to evaluate the significance of the trend of field 
strength with the phase of the solar cycle (-94Gyr _1 ). 

First experiment: Assume that the maximum magnetic field 
strength of sunspots has a constant distribution at any moment of 
a solar cycle, i.e., the PDF of the maximum field strength, B max , 
is the same for each year and is given by Fig. [7j The number 
of sunspots close to a solar maximum is about 100 times more 
than at the solar minimum (cf. Fig. [8j. Therefore, the statistics 
of sunspots around the time of maximum activity is better and 
will sample the PDF correctly, whereas close to minimum activ- 
ity strong deviations from the true PDF may occur because of 
the sparse sample. This can lead to an apparent temporal trend 
in field strength caused purely by the variation of the sample 



size, where both a decrease or an increase in the field strength 
can result. To reject the hypothesis that the measured system- 
atic variation of -94 G yr~' is caused only by the large number of 
sunspots close to the maximum of a cycle and the low number 
close to the minimum (Fig. [8]), we performed the following test. 
Figure [8] shows the yearly sunspot number from the Solar 
Influences Data Center (SIDC) during cyc le no. 23 and the early 
years of cycle no. 24 dSIDC-teamll201 lb . For each year from 
1999 through 2009, we created a sample of B max , according to 
the observed distribution of field strengths (Fig. |7). The size of 
the sample population was set to the annual sunspot number, 
making the PDFs of the the field strength in individual years 
roughly identical within the statistical fluctuations. We then an- 
alyzed this synthetic temporal sample in the same manner as 
the real data. Using both the Bayesian (Sect. l4~2l and the least- 
square methods, we fit a straight line to measure any change of 
the mean of B max . The test was repeated 100 times. The average 
and rms of the resulting slopes were zero and 1 1 G yr" , respec- 
tively. The two methods resulted in identical average and rms 
slopes. The sample of 100 slopes had a maximum and minimum 
of +23 and -26 G yr~ ' , respectively. We note that, if one plots all 
such artificial solar cycle curves, the width of the region covering 
the 100 curves increases toward the end of the cycle as expected 
because of the poorer statistics. That means the upper bound- 
ary for statistical fluctuations of less than 30Gyr~' as the sole 
reason for the observed systematic trend is significantly smaller 
than the slope we have derived from the observations. 



Second experiment: Although the first test showed that the 
trend is real, it is not clear how far it is influenced by the statis- 
tical fluctuations. We performed a second numerical experiment 
to elaborate on this issue. Like the first test, we created artifi- 
cial samples of B max . The sample size for each year was as in 
the first test. The PDF of the B max was, however, different for 
each year in the sense that the mean value of the PDF was de- 
creased by -94 G yr -1 . Because of the lower number of sunspots 
at the end of the cycle, it is then possible to get a slope that is 
too flat or too steep compared to the actually imposed slope. The 
number of repetitions and the analysis method were identical to 
the first experiment. The average and rms of the slopes derived 
in the second experiment were -94 and 12 Gyr -1 , respectively, 
both for the Bayesian and least-square methods. The maximum 
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and minimum slopes of the sample were -68 and -130 Gyr , re- 
spectively. That means that with a complete sample the measured 
slope can be off by about 36 G yr" 1 in a worst case scenario. 



Third experiment: The second experiment showed that, if there 
is a real trend in B max with the phase of the solar cycle, it is mea- 
surable despite the different statistics of sunspots at the maxi- 
mum and minimum of the cycle. However, we do not have all of 
the sunspots in our sample, only those that were observed with 
TIP. This can introduce a significantly larger uncertainty in the 
derived slope compared to the second experiment. Therefore, we 
repeated the second numerical test with one change: instead of 
using the annual sunspot number of SIDC, we used the statistics 
of the TIP sample. The resulting slopes then spanned a range be- 
tween -130 and -66 G yr~ : with an average of -84 and a standard 
deviation of 9 Gyr -1 . This means the statistics of the TIP sam- 
ple has about the same amount of uncertainty as the complete 
sample with a caveat that it underestimates the correct trend by 
lOGyr -1 , well within the \cr range of 14Gyr~'. This ensures 
us that the slope of -94 G yr~' is significant and measurable with 
the statistics of our sample. 



5. Discussion 

5.1. Tracers of the solar cycle 

The number of sunspots and sunspot groups ar e the most com- 
mon indicators for studying the solar cycle (lHathawav et al.l 
l2002l and references therein). Nevertheless, other parameters 
with more straightforward physical meanings show similar, if 
not better traces of cyclic variations. These parameters, such as 
the magnetic field strength of sunspots, usually have a caveat: 
their systematic recording started only in the early 20th century. 
In other words, they cannot compete with the 400-year record 
of sunspot numbers. Recent investigations of the sunspot num- 
bers also indicate that large sunspots correlate better with chro- 
mospheric and cor onal indicators of the solar cycle such as the 
10.7 cm radio flux (iKilciket alJl201lh . 



5.2. Continuum intensity vs. magnetic field strength 

Our results reproduce earlier observations regarding the 
general behavior of the maximum magnetic field and mini- 
mum continuum intensity of sunspots: the greater the field 
strength, the lower the continuum inte nsity. Our Fig . |3 



compares reasonably well with Fig. 2 of Livingston! (1200: 
Similar curves of the relation between intensity and field 
strength for this and other wavelength ranges a r e given in 
literat ur e (e.g., Martinez Pillet & Vazaue A 119931: Lites et al. 
19931: iKeppens & Martinez Pilletl 1 1 996k IStanchfield et al. 



1997t iNorton & Gilmanl 120041 : Leonard & Choudharvl I2008L 
Wesolowski et al. 2008; Schad & Penn 2010), where those of 
Gurman & House (1981) and Livin gston! (|2002|) match ours 
closest because of the pronounced linear behavior. In some 
of the works the relation between intensity and field strength 
was derived inside a single sunspot, not for a sunspot sample, 
sometimes also yielding clearly nonlinear relations. 

5.3. Temporal evolution of 6 max 

Solar cycle no. 23 was weake r than cycle no. 22 (e.g., 
lAtac & Ozgu^l 120061) . Livingston! d2002l) report that sunspots 
in cycle no. 23 had an excess of smaller and brighter um- 



brae compared to cycle n o. 22. Using 1.56// m Kitt Peak data, 
iPenn & Livingston! (1201 lh suggest that the strengths of sunspots 
has steadily decreased since the last solar maximum. Although 
our finding (Fig.|6]l in the declining phase of cycle no. 23 is con- 
sistent with them, we do find an increase in B max since the be- 
ginning of the solar cycle no. 24. As seen in the lefthand panel 
of Fig. [6] the average field strength of the years 2010 and later 
are well above the linear trend that indicates the extrapolation of 
the decrease during the declining phase of c ycle n o. 23. That is 
in contrast to Fig. 1 of IPenn & Livingston! (1201 ll) . who predict 
a continuous decrease in the sunspot field strength and conse- 
quently on the solar magnetic activity across cycles. We note, 
however, that the field strength at the start of cycle no. 24 seems 
to be slightly lower than at the same stage in cycle no. 23, and the 
continuum intensity correspondingly slightly higher. This will 
have to be confirmed by additional data in the following years. 

The weakening rate of the magnetic field strength in the 
descending phas e of cycle no. 23 i n the TIP data is about - 
94(± 141Gvr~'. |Penn & Livingston! (120061) reported SlGyT 1 
The size o f our sample is smaller th an that of Livingston! (|2002|) . 
as well as lPenn & Livingston! (120061) . However, the precision of 
the measurements, and the calibration accuracy is superior to the 
mentioned references because of using vector-polarimetric data 
instead of intensity profiles. 

IWatson et al.l (1201 ll) used Soho magnetograms and analyzed 
more than 30,000 sunspots. They reported a systematic decrease 
of -70 G yr in the declining phase of the cycle no. 23, in agree- 
ment with our results. These results should be taken with care 
because the Soho/MDI instrument cannot accurately measure 
the magnetic field in sunspots because of lin e saturation effects 
dBerger & Litesll2003l) . iPevtsov et ail (1201 lb compiled observa- 
tions of seven (former Soviet) solar observatories. They report a 
rate of -119 Gyr _1 for cycle no. 23, which is the strongest rate 
in the four cycles they studied. Despite some (systematic) un- 
certainty in part of their data around 1965, they did not find any 
significant secular trend in field strength over four and a half so- 
lar cycles. 



5.4. Temporal evolution of umbral intensity 



lAlbregtsen & Maltbvl (1 19781) and lAlbregtsen et all (1 19841) found 
that the amount of variatio n in the umbra l intensity during a cy- 
cle is about 10% (see also. iMaltbvll 19921). Our resu l ts are there- 
fore consistent with their analysis. iMathew et al.l (120071) used 
Soho data and did not find any significant enhancement in the 
umb ral intensity as a function o f the solar cycle phase. In con- 
trast, |Pennj&MacDcmald| (|2003) find a cycle-dependent oscilla- 
tion of the umbral intensity. Our data show an increase in the 
umbral continuum intensity with the phase of the solar cycle. 
We find a slope of 1.3 (+0.4)% I c (1.56^m)%yr 1 in the TIP 
data of cy cle no. 23. Our results are consistent with the obser- 
vations of lAlbregtsen et all d!984l) . who reported an increase of 
1.2 and 1.9 % per year in cycles no. 20 and 21, as well as with 
IPenn & Living ston ( 2006) who find an increase of 1 .9 % per year 
from 2000 to 2006. 



5.5. Umbral area 



iBrants & Zwaanl (1 19821) found t hat larger umb r ae ha ve a 
high er field strength, sim i lar to iKopp & Rabin! (1 19921) and 
iMartinez Pillet & Vazquez! (Il993l). We find a sim ilar trend in 
our results (Fig. [5). Recently, IWatson et al.l (1201 ll) report a de- 
pendence of the area of the largest sunspots on the solar cy- 
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cle. This, however, doe s not imply any trend in th e size of in- 
dividual sunspots. Like iPenn & MacDonaldl (120071) . we did not 
find any systematic trend between the phase of the solar cycle 
and the sunspot size. These observations therefore hint at a con- 
stant distribution of the s unspot area, e.g., a log-normal distri- 
bution with a fixed slope dBogdan et alj|1988l) . The larger num- 
ber of sunspots close to the maximum of a cycle naturally in- 
creases the number of large sunspot groups. Because these large 
sunspots are located at the tail of the distribution, they are prone 
to larger statistical fluctuations compared to other quantities re- 
sulting from all sunspots such as B max . The lack of a system- 
atic change of sunspot area is still puzzling because we mea- 
sure a difference of some 700 G between the average B max in 
the maximum and minimum of cycle no. 23, comparable to the 
600 G reported by IPenn & Livingston! (120061) . Taking Fig. |5] at 
face value shows the relation between field strength and area, 
a rate of -28arcsec 2 yr~' would be expected. A least-square fit 
to the average umbral areas per year resulted, however, only 
in a slope of 2(±4)arc sec 2 yr~ 1 . It was already mentioned by 
iBrants & Zwaanl dl982l) th at the scatter between B max and um- 
bral size is partly intrinsic. IPenn & MacDonald (12007b suggest 
that the relation between the B max and umbral area might vary 
during a solar cycle. Therefore, it remains an open question to 
verify whether the diversity of physical conditions in different 
umbrae vary with the solar cycle in a sense that it prevents a 
cyc lic trend of umbral si ze. 

iMathew et al.l d2007l) report a clear trend between the um- 
bral intensity and the umbral size in Soho data (see also 
IWesolowski et al.l|20 08). Although the umbral contrast is lower 
at 1.56/im, there is a similar trend in our data as well. As 
seen in Fig. |4] larger sunspots have a lower continuum inten- 
sity than smaller sunspots, in agreement wi th previous studies 
( IBrants & Zwaarll982tlKopp & Rabinlll992h . 



5.6. Implications for the solar irradiance 

An excess or lack of larg e dark umbrae changes the total so- 
lar irradiance (TSI, e.g., ISteinegger et al.l fl990b . Our finding 
thus also has some implications for the TSI variations in solar 
cycle no. 24. There are l ong term variations i n the solar ac- 
tivity and irradianc e (e.g., iTapping et al.l 120071: iFrohlichl 12009: 
IShapiro et al.ll2.01 lb . In particular, because the Sun apparently 
had an unusually active period in the last couple of decades 
(lUsoskin et alJl2003c iBonev et al.ll2004l) . we cannot completely 
rule out the existence of a long term trend. However, in spite 
of suggestions that cycle no. 24 and the pre ceding activity mini- 
mum are very different from previous on es (iTripathv et al .120101 
iJain et al.l20 lit (Tapping & Valdesl201 lb . the influence of a long 
term trend seems to be less than the cycle-dependent pattern on 
the timescale of a solar cycle. 

6. Conclusion 

A sample of sunspot maps observed with the Tenerife Infrared 
Polarimeter was analyzed to elaborate on a possible dependence 
of the magnetic field strength of sunspots on the phase of the 
solar cycle. The vector-polarimetric data enabled us to use full 
Stokes measurements to estimate the variation in B max rather 
than only using intensity profiles as used in case of Kitt Peak 
data. We found a significant decrease in the measured mag- 
netic field strength from the maximum to the minimum of the 
sunspot cycle no. 23 with a rate of -94 (+ 14) Gyr -1 . At the 
same time, the minimum umbral continuum intensity increases 
by 1.3 (+ 0.4)% of I c yr" 1 , while there is no significant trend of 



the umbral area with the phase of the solar cycle. The magnetic 
field strength of sunspots increased at the beginning of solar 
cycle no. 24. This indicates that the observed trend of B max is 
likely to be a cyclic pattern rather than a long-term evolution. 
Comparing the field strength in the early stages of cycle nos. 23 
and 24, we find a slight reduction in B max in the new cycle that, 
however, will need to be confirmed by additional data in the fu- 
ture. 
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